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Abstract 

We present a new polar representation of quaternions inspired by the Cayley-Dickson represen- 
tation. In this new polar representation, a quaternion is represented by a pair of complex numbers 
as in the Cayley-Dickson form, but here these two complex numbers are a complex 'modulus' and a 
complex 'argument'. As in the Cayley-Dickson form, the two complex numbers are in the same com- 
plex plane (using the same complex root of —1), but the complex phase is multiplied by a different 
complex root of —1 in the exponential function. We show how to calculate the amplitude and phase 
from an arbitrary quaternion in Cartesian form. 

1 Introduction 

It is well-known that the complex exponential e 10 — cos 9 + % sin generalises to quaternions by replacing 

1 by any unit pure quaternion fi, since any unit pure quaternion is a root of —1 [G, pp. 203, 209] [ , §167, 
pl79]. Hence, any quaternion may be represented in the polar form q — |g|e M6> where 8 is a real angle [8, 
§2.3] [I, §12.7]. A small difference between the complex and quaternion cases is that in the quaternion 
case the argument 9 is conventionally confined to the interval [0, ir}. This is because the modulus of the 
vector part of the quaternion is always taken to be positive (there is no convenient way to define an 
orientation in 3-space which would permit the sign of the vector part to be determined). 

The classic polar form just discussed is not the only possibility for representing a quaternion. There 
are representations based on Euler angles using three exponentials in which the arguments of the expo- 
nentials are directly linked to the Euler angles. Billow, for example, in [2, Theorem 2.3] and [3, §V.A] 
quotes the following formula, and gives an algorithm for finding the three angles: 

q = \q\e^e k ^e je 

See also [1, §12.11] for a discussion of Euler angles and quaternions. 

The polar form we present in this paper is different to both of the preceding forms, in that it represents 
a quaternion using a single exponential with a complex argument and a complex 'modulus'. 

2 Polar form 

We now show that every quaternion admits the following polar form, in addition to those presented in 
the previous section: 

q = Ae Bj 

where A — a + bi and B = c + di are complex and a, b, c and d are real. Before proceeding, we note 
the analogy between this polar form, based on two complex numbers, and the Cayley-Dickson form of a 
quaternion q = (w+xi) + (y+zi)j — w + xi+yj + zk which is also based on two complex numbers. In the 
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Cayley-Dickson form the first complex number has real and imaginary parts which are simply the first 
two parts of the quaternion in Cartesian form, and the second complex number has real and imaginary 
parts which are simply the third and fourth parts of the quaternion in Cartesian form. There is a subtlety 
to the Cayley-Dickson form which we must mention before proceeding, otherwise what follows may not 
be clear. The subtlety lies in the difference between a complex number such as a + bi and a quaternion 
a + bi + Oj + Ofe. In this paper, as in the classic Cayley-Dickson form, the complex number is treated 
as a degenerate quaternion and we handle complex numbers algebraically as quaternions in which two 
components happen to be zero. The usual rules of quaternion algebra therefore apply, including the rules 
i 2 = j 2 = k 2 — ijk = —1. 

The polar form presented in this paper is very subtle. It shares with the Cayley-Dickson form the 
idea of a construction based on two complex numbers in the same complex plane with % as the root of 
— 1, the second of these (the argument of the exponential) being multiplied on the right by j, just as in 
the Cayley-Dickson form. Unlike the Cayley-Dickson form, the parameters of the polar form presented 
here are not trivially obtainable from cither the Cartesian form of a quaternion or the classic polar form. 

In what follows, we first need to prove that e Sj is a special type of quaternion, with only three 
components: a + 7ji + 5k, that is, it has a zero coefficient of i. 

Lemma 1. Given an arbitrary quaternion in the form p = (c + di)j = cj + dk, its exponential is given 
by: 

p 

e p = cos \p\ + — sin \p = a + jj + Sk 
\P\ 

Proof. Represent p by the product of its modulus and a unit quaternion obtained by dividing by the 
modulus: 

Comparison with the classic polar form of a quaternion given in §1, namely e^ e — cos# + ^tsin#, where 
fj, is a unit pure quaternion, is sufficient to prove the result, since p/\p\ is a unit pure quaternion which 
may be identified with fi, and \p\ is real and may be identified with 6. Thus a — cos \p\, 7 = (c/|p|) sin \p\ 
and S = (d/\p\) sin \p\. □ 

Now we present a polar form as previously outlined with a complex 'modulus' and a complex argu- 
ment, and we show how to express an arbitrary quaternion in this form, that is, how to find the complex 
'modulus' and argument. 

Theorem 1. Every quaternion q = w + xi + yj + zk, where w, x,y, z f= R can be expressed in the form 
q = Ae B i , where A,B £ C, specifically A = a + bi and B = c+ di. 

Proof. The proof is by demonstrating that we can find A and B without any constraints on q. Clearly, 
if q = 0, A must be zero, and the value of B is unimportant. Otherwise, we note that if we put the 
complex 'modulus' A in polar form as A — \A\e ei , we can immediately see that \A\ = \q\. Therefore 
in what follows we may assume, without loss of generality, that we are working with quaternions q of 
unit modulus, and that therefore A = a + bi is of unit modulus. Therefore we have, using the result in 
Lemma 1: 

q = Ae 3 -* = (a + bi)(a + jj + 8k) = aa + bai + (07 — bS)j + (aS + lr/)k 

Now it is possible to see that we can obtain a + bi quite trivially, provided a ^ 0. Assuming a ^ 
for the moment, and since a + bi has unit modulus, all we need is to take a complex number formed 
from aa and ba (the first and second components of q, which are the scalar part and the coefficient of 
i), and normalise this complex number by dividing by its modulus. However, there is a small problem: 
there is an ambiguity of sign. We do not know the sign of a, and this means we have not determined 
a+bi completely since — (a + bi) is equally valid. It is clear that this ambiguity is insurmountable, since 
negating a + bi can be compensated for by negating a + 7J + Sk. 

In the case where a — 0, but \q\ > 0, the first and second components of q will be zero. Thus, q is of 
the form 7.7 + Sk and we may take a = 1 and 6 = (thus A — I). The factorization here is not unique, 
and in §3 we present an example to show this (we could have chosen a = and 6=1, and thus A = i). 
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Once we have obtained a + bi it is simple to find a + jj + Sk by dividing g on the left by a + bi, or 
equivalently, by multiplying on the left by the inverse 1 of a + bi. 
To obtain Bj we can use the quaternion log function: 

Bj = log (a + 7ji + Sk) 

Numerically this is a simple process, since the quaternion logarithm can be computed using a complex 
logarithm operating on a complex number isomorphic to the quaternion. The resulting complex value 
can then be used trivially to construct a quaternion with vector part oriented in the same direction as 
the original quaternion 2 . 

It should be clear from Lemma 1 that Bj must be of the form cj + dk and therefore the real and 
imaginary parts of B are trivially obtained from the second and third imaginary components of Bj (the 
coefficients of j and fe) . □ 



2.1 Sign ambiguity 

We noted above that the sign of A has an ambiguity which can be compensated by an alteration to B 
but we did not discuss what form this alteration would take. We merely pointed out that negating A 
required a + 77 + Sk to be negated. 

It is easy to see that negation of a + -fj + Sk may be accomplished by adding ir to the modulus of 
B. From Lemma 1 adding n to \B\ means that the angles in the cosine and sine are augmented by tt 
which negates their results. It is not so easy to see that the converse is not the only possibility. In fact 
alternative solutions exist that modify B to negate a + jj + Sk, and the algorithm presented above for 
finding B using a quaternion logarithm will yield these solutions. For the present paper we confine our 
discussion of this issue to some numerical examples. 



3 Numerical examples 

Let q = 4(1 + i+j + fe). Then |<7| = 1 and aa = ba — ~ which gives A = ^= + -^i. Multiplying q on 
the left by the conjugate of A gives + A^j and the quaternion logarithm gives Bj = irj /4 (therefore 
B is real in this instance). To verify the result: 

/I 1 7T . , 1 1 ... 7T . . 7T. . 1 1 ... 1 .1,1,,..., 

(—7= H 1=1) exp —3 — ( — H — i)(cos - + j sin -) = (— H —i)( — +3—=) = + i + 3 +k) 

For a second example, we take q = 1 + 2i + 3j + 4fc. In this case we state the result obtained 
numerically which is A = 2.4495 + 4.8990i and B = 1.1317 - 0.2058i. The modulus of B = 1.1503. 
Hence (apart from rounding errors caused by the limited precision of the intermediate results shown 
here) : 

(2.4495 + 4.8990i) exp(1.1317j - 0.2058fe) = 
(2.4495 + 4.8990i) (cos 1.1503 + (0.9839j - 0.1789fc) sin 1.1503) = 
(2.4495 + 4.8990i) (0.40825 + 0.89815.? - 0.1633A;) = 1 + 2i + 3j + Ak 

Finally, we show an example where the value a in Theorem 1 is zero, and the factorization is not 
unique. Consider the factorization of the quaternion with value k. The following two factorizations will 
yield this value. Firstly A — i and B = ^: 

i exp (^-j) = i (cos | + j sin |j = i(0 + j) = k 
Secondly, A = 1 and B = ^i: 

/7T..\ 7T 7T . . 

1 exp ^-13 J = cos - + %3 sin - = + 13 = k 

1 The inverse of a + bi is a — bi if a + bi has unit modulus (as assumed) since the inverse is the conjugate divided by the 
modulus. 

2 See, for example, the file log.m in our Matlab® quaternion toolbox [7] 
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